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Abstract 

Phase transitions can occur in one-dimensional classical statistical mechanics at non-zero 
temperature when the number of components N of the spin is infinite. We show how to solve 
such magnets in one dimension for any N, and how the phase transition develops at N = oo. 
We discuss SU(N) and Sp(N) magnets, where the transition is second-order. In the new 
high-temperature phase, the correlation length is zero. We also show that for the SU(N) 
magnet on exactly three sites with periodic boundary conditions, the transition becomes first 
order. 

1 Introduction 

It has long been known that phase transitions are uncommon in one-dimensional classical statis- 
tical mechanics. An old argument by Peierls shows that in models at non-zero temperature with 
local interactions and a finite number of degrees of freedom, order is not possible: the entropy 
gain from disordering the system will always dominate the energy loss. There are (at least) three 
ways of avoiding this argument. The first two are well understood. A system at zero temperature 
can of course order: the system just sits in its ground state. A system with long-range interac- 
tions can have an energy large enough to dominate the entropy. In this paper, we will discuss 
in depth a third way of obtaining a phase transition in one dimension. This is to study systems 
with an infinite number of degrees of freedom per site. 

In particular, we will study magnets with O(N), SU(N) and Sp(N) symmetry. We will 
see that there can be a phase transition in the N — > oo limit. We solve these one-dimensional 
classical systems for any N, and show how the transition occurs only in this limit; for finite N 
all quantities depend on the temperature analytically. The infinite number of degrees of freedom 
has roughly the same effect of increasing the effective dimensionality, but the phase transition is 
very different from those in higher dimension. It is not a phase transition between an ordered 
phase and a disordered one, but rather between a disordered phase and a seriously-disordered 
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one. In the seriously-disordered phase, the system behaves as if it were at infinite temperature. 
The entropy has dominated the energy to the point where the energy term does not affect the 
physics; each spin is effectively independent. The infinite number of degrees of freedom means 
that this serious disorder is possible even at finite temperature. 

The paper is a companion to one by Tchernyshyov and Sondhi [j]]. There it is shown that 
in some magnets, a mean-field calculation yields a phase transition in any dimension. Since 
mean-field results are exact at N — > oo, this predicts the phase transition we observe here. Their 
computation also predicts that there is a first-order phase transition for the SU(N — * oo) magnet 
on just three sites with periodic boundary conditions. Remarkably, this first-order transition 
happens only for precisely three sites; for any other number of sites greater than 1 there is a 
second-order transition. 

It has long been known that phase transitions can occur as N — » oo in zero-dimensional 
matrix models P] . Phase transitions in one dimension at infinite N were studied in M, |J] . In 
particular, the largest eigenvalue for the SU(N) case discussed here was computed in || for any 
N. Here will develop the necessary techniques systematically, and extend these results in several 
ways. We explicitly find all the eigenvalues of the transfer matrix for these magnets. All these 
results are completely analytic in N and in the inverse temperature (3 as long as N is finite. The 
singularity and a phase transition can develop when N — > oo and (3 — > oo with (5/N remaining 
finite. Knowing all the eigenvalues and their multiplicities explicitly for any N lets us show that 
there can be a phase transition as N — * oo even for a finite number of sites in one dimension. 

In section 2, we find all the eigenvalues (and their multiplicities) of the transfer matrices 
in a variety of one-dimensional magnets. In section 3, we use these results to study the phase 
transitions which occur as the number of sites and N go to infinity. Most of these phase transitions 
are ferromagnetic, but one is antiferromagnetic. In section 4, we study the first-order transition 
for the three-site SU(N) chain. In an appendix we collect some useful mathematical results. 

2 Solving one-dimensional magnets at any N 
2.1 The rotor 

To illustrate the procedure, we start with a simple rotor, the classical XY model in one dimension. 
The spin is defined by a periodic variable < 9 < 2tt, and spins 9 and 9' on adjacent sites have 
energy 

E o(2 )(9,0') = -cos(9-9>). (1) 
To compute the partition function of this system, define a transfer "matrix" 

T(9,9') = e-^ e ' e '\ 

Since the variables of the system take continuous values, this isn't really a matrix, but rather the 
kernel of an integral operator. It takes functions of 9 to functions of 9' by 

9(0') = J o ^T(9,9')f(9'). 

To compute the partition function, we need eigenvalues Xi of T. Because the spins take values 
on a compact space (the circle here), the eigenvalues are discrete and hence labeled by a discrete 
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index I. The corresponding eigenfunctions fi{9) obey 



kfi(o') = j Q ^ne^e'W). (2) 

For the energy (|l|), the fi(9) are obviously 

fl = e* W . 

The index I must be an integer to preserve the periodicity under 9^9 + 2tt. To see that these 
are eigenfunctions, note that 

\ U6 _ ^_ 0cos(6-6') il(e') 

2n 











27r rift' 

au e /3cos(6»') e ii{e~e') 
2n 



The integral then can be evaluated for any I in terms of a Bessel function: 







= IM- (3) 
The partition function for L sites with periodic boundary conditions is then 



When L is large enough, the sum is dominated by largest eigenvalue, which here is the I = 
state. The internal energy density U of the system is then 

u °^--idi3 lnZ --iWy 

All other quantities such as correlators can easily be found as well, since we have an explicit and 
complete set of eigenvalues and their multiplicities. 



2.2 O(N) magnets 

The eigenvalues of the 0(2) problem are found by Fourier transforming the transfer matrix. 
What we need to do for more general cases can be summarized as Fourier analysis on manifolds 
more general than the circle. In other words, we want to expand a function taking values on a 
manifold Ai into a series, e.g. 

T(M)=J2\ l fi(M) 

i 

where the fi(M) are complete set of orthonormal functions. The eigenvalues of the transfer 
matrix are the coefficients of the expansion in this basis. 

The problem of Fourier analysis on all the manifolds of interest has been solved already. The 
key is to exploit the symmetry. A familiar example is where the spins take values on the two-sphere 
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S 2 , where the eigenfunctions are called spherical harmonics. In coset language, the two-sphere 
can be described as the manifold S 2 = 0(3)/0(2): the 0(3) group consists of rotations, while 
the 0(2) subgroup is the set of rotations which leave a given point invariant. Thus different 
points on S 2 take values in 0(3)/0(2). We parameterize the two-sphere by the usual spherical 
coordinates: a unit three- vector s = (sin ^ sin ^, sin#cos</>, cos 9). To make progress, it is crucial 
to consider an energy invariant under the 0(3) rotation group, namely 

Eq(3) = —s- s' = — cos#cos#' — sin 9 sin 9' cos (<p — <f>'). 

We can expand the transfer matrix energy into irreducible representations of the rotation group, 
labeled by an angular momentum I and an L z component m. The eigenfunctions of the transfer 
matrix are expressed in terms of Legendre polynomials Pi, whose explicit definition will be given 
below. One way of showing that the Legendre polynomials are eigenfunctions of the transfer 
matrix is to show that they obey an addition theorem, for example 

Pi(cos8 cos6' + sin9 sm9' cos(cj) — <fi')) 

= P ; (cos 9)Pi(co S 9') + 2 V f ~ m |; ^ m (cos 9)P lm (cos 9') cos(0 - <f>'), 

where the P/ m are called associated Legendre polynomials, with Pi = Piq. One can then expand 
the function e^ x in terms of Legendre polynomials, and then use the addition theorem to split 
the 9 and 9' dependence. This leaves an integral for the eigenvalue. 

To just obtain the eigenvalues of the transfer matrix, one does not have to go to all the 
complications of generalized addition theorems. Mathematicians have developed more in-depth 
ways of deriving the eigenfunctions, and then the addition theorem comes as a consequence of 
the computation. A geometric method is discussed in []|], while a much more explicit method is 
discussed in ||. We will require the methods of the latter in order to treat the Sp(N) magnet, 
where the spins do not take values on a symmetric space. (A symmetric space G/H has H a 
maximal subgroup of G; the importance here is that when the spins take values in a symmetric 
space, the transfer matrix depends on only one parameter.) 

First we find the eigenvalues for an 0(iV)-invariant magnet, where the spins s take values on 
the manifold 0(N)/0(N - 1), which is the (N - l)-sphere S N ~ 1 . We take the energy between 
nearest neighbors to be 

E {N) = -s- s>. (4) 
An eigenvalue A; associated with eigenfunction fi(s) is given by the equation 

>H iN) fm = I [Ds']e^ sl fi(s'). (5) 
Js 1 *- 1 

where [Ds'] is the usual measure on the N — 1 sphere, normalized so that J[Ds'] = 1. The unit 
iV- vector s depends on JV- 1 angles, but the energy only depends on the angle between the two 
spins. This means that to compute the eigenvalue, we need do only one integral. Put in terms of 
the 0(3) spherical harmonics, it means that the eigenvalues depend only the value of I and not 
m. Explicitly, if one sets s = (0,0, ... ,0, 1), the integrand in (|5|) depends on only a single angle 
9. We can then do the integral over all the other angles in (|5|), and the measure reduces to Q 
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The eigenfunctions for all the magnets we study here can be written in terms of Jacobi 
polynomials P, (x). These are orthogonal polynomials of order x l , and a number of useful 
properties are collected in the appendix. The eigenfunctions for the O(N) magnet with energy 
(141) are given by P( r ' r )(cos60, where r = (N - 3)/2. These are often called Gegenbauer or 
ultraspherical polynomials. They are indeed orthogonal with respect to the measure The 
eigenvalues are then given by 

1 p (1) "r((iv-i)/2)r(i/2) J desm {e)e p {cos9) 

Using the integral in the appendix for \x = v= (N — 3) /2 gives 



,0(N) _ „_ T(N/2)T(N/2-l/2 + l) 

r(i/2)r(jv-i + 2i) 



x yW) = e -p w Z-A_L — . / ^ J M(N/2 - 1/2 + I, N - 1 + 1, 2/3). 



The function M is called Kummer's function, and is a confluent hyper geometric function. Its 
definition and the differential equation it satisfies are given in the appendix. Using the double- 
argument formula for gamma functions || and the identity p9| ) with v = N/2 — 1 + 1 gives our 
result 

/ a\ N/2-1+1 

X ? (N) = T(iV/2) In/2-iM (7) 

Note that the eigenvalues reduce to the rotor result (||) when N = 2. Note also that O(N) 
ferromagnets and antiferromagnets are essentially the same, because redefining s — * — s for every 
other spin sends (3 — > —(3. This transformation 6 — > ir — 9 leaves the measure of the integral 
invariant, and \? {N) (-(3) = (-l) z Af (iV) (/3). 

As a function of /3, the Bessel function is analytic for all (3. Moreover 1^(0) > Xy+i(/3) for 
any f3 as long as 7 is positive. Thus the free energy does not have any singularities as long as 7 
remains finite. One might hope something interesting happens when N — ► 00, but in Section || 
we will show that in this case there is still no transition. 



2.3 SU(N) magnets 

The computation for the SU(N) case is very similar to that of the 0(N) case, but we will see in 
the next section how there is a completely new phase. 

The SU(N) magnet is defined in terms of a complex A^-vector z obeying z* ■ z = 1. The 
energy for adjacent sites is 

E su W=-\z*-z'\ 2 . (8) 

This energy is not only invariant under global SU(N) rotations, but is invariant under local 
(gauge) transformations z{x) — > e ta<yX ^z{x) at any site x. The vector z takes values on a complex 
sphere U(N)/U(N — 1), which as a manifold is identical to the real 2A" — 1-sphere 0(2N)/0(2N — 
1). However, the U(l) gauge symmetry can be used to effectively reduce the number of degrees 
of freedom in the problem by 1. For example, one can set the last component of z to be real at 
every site, so z effectively takes values on the manifold 

U(N — 1) x U(l)' 

This manifold is a symmetric space, and is known as the complex projective space CP N ~ l . 
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Like the O(N) case, the energy between adjacent sites depends only on a single angle 8, where 
z* ■ z' = cos 9 (any phase can always be gauged away). Then the normalized measure can be 
written as [||, |j 

f [Dz'\ = 2{N - 1) f d9 sin 2Af - 3 9 cos 9 (9) 

Jap"- 1 JO 
The eigenvectors of the CP N_1 transfer matrix are Jacobi polynomials as well, namely || || 

The N = 2 case reduces to the ordinary real two-sphere, and indeed the P^ ' ' are ordinary 
spherical harmonics. The SU(N) symmetry means that the eigenvalues have degeneracies. For 
example, we saw in section 2.1 that for N = 2 the eigenvalue depends only on I and not the S z 
value m. For N = 2, the degeneracy Di of Xi is of course 21 + 1; the generalization to all N is || 

= (2/ + jv-i)r 2 (/ + jv-i) 

' r(iv)r(iv-i)r 2 (/ + i) 1 ' 

The eigenvalues are given by 

»7r/2 



x SU(N) p (N-2,o) {1) = 2{N _ 1) r d e S m 2N - 3 eco S 9e^ os2e p[ N - 2fi \co S 2e) 

Jo 

dx(l - x) N - 2 eft 1+ *V 2 Pl N - 2 >°\x) (11) 



JV-1 



2 N-1 

where x = cos 29. Using the integral in the appendix for fi = N — 2, z^ = gives 

x su W = ^ rgvyp + i) M{1 + % N + ^ (12) 

For finite iV and /?, these are analytic functions. Note that, unlike for O(N) magnets, (3 > and 
/3 < are not equivalent here. For positive (3, we have a ferromagnet: the energy favors aligned 
spins: a neighbor of z = (0, 0, . . . , 0, 1) is in the same state z' = z (up to a phase). However, 
(3 < does not correspond to an antiferromagnet: the energy favors the neighbor z 1 being in 
any state orthogonal to z, i.e., z* ■ z' = 0. There is a unique orthogonal state only for N = 2. 
Thus, with the exception of SU(2), the SU(N) model (||) is ferromagnetic when /3 > 0, but is 
not an antiferromagnet when (3 < 0. For a proper generalization of an antiferromaget, one needs 
to designate among the N — 1 vectors orthogonal to z = (0, 0, . . . , 0, 1) one that can be called 
"antiparallel" to it, for example, z' = (0, 0, . . . , 1, 0). Doing so manifestly breaks the SU(N) 
symmetry, as we discuss next. 

2.4 Sp(N) magnets 

The SU(2N) magnet can be deformed in an interesting way by breaking the SU(2N) symmetry 
to Sp(N). The main novelty of this Sp(N) magnet is that the phase transition can occur at an 
antiferromagnetic value of (3. The energy between adjacent sites is 

=- a \z* -z'\ 2 -b\zJz'\ 2 . (13) 

where z is a 2iV-dimensional complex vector, and J is the 2N x 2N dimensional matrix ia y ® I, 
where I is the iV-dimensional identity matrix, and a x the Pauli matrix. The second term breaks 
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the SU(2N) symmetry, but preserves an Sp(N) subgroup. The case b = reduces to the SU(2N) 
ferromagnet discussed above. The case a = 0, studied by Read and Sachdev [f?J, is a large- N 
generalization of an antiferromagnet. Indeed, the energy 



\zJz'\ 2 



N 



t=l 



is minimized when vectors z and z' are locked in adjacent flavors 2i — 1 and 2i, e.g., z = 
(0, 0, . . . , 0, 1) and z' = (0, 0, ... , 1, 0), analogues of spin-down (0, 1) and spin-up (1, 0) in SU (2). 

Like the SU(N) magnet, the generalized model has a local U(l) symmetry. Just like the 
complex sphere U(N)/U(N — 1) is identical as a manifold to the real sphere 0(2N)/0(2N — 1), 
the "quaternionic sphere" Sp(N)/ Sp(N — 1) is identical as a manifold to the complex sphere 
U{2N)/U(2N -I) (which in turn is equivalent the real 47V-l-sphere 0(4N)/0(4N- 1)). When 
we fix a gauge, then z takes values on the manifold 

Sp W (14) 



Mi 



Sp{N - 1) x 17(1) 

It is easiest to study first the special point a = b, where the gauge symmetry is enhanced to 
SU(2) = Sp(l). The SU(2) gauge symmetry (a subgroup of the Sp(2N)) mixes the 2i — 1 and 
2i components of z and also mixes z and z* . Precisely, if we arrange these components into the 
matrix 

' 22i-l 

z 2i *5i_i 

then all the Mj transform under the gauge symmetry as 

Mi -» MiUi, 

where Ui is an element of SU{2). The energy is invariant under these transformations even though 
Ui can be different at every point. If we use this symmetry to fix a gauge (e.g. zn-i = sm(0jv-i)> 
zn = cos(^at_i)), then z for a = b takes values on the manifold 

Sp(N) 



Sp(N - 1) x Sp(l) ' 



which we call QP 1 *' 1 . Solving this case is almost the same as the SU(N) case, because this is 
a symmetric space. The energy between adjacent sites depends on only one variable, and the 
measure becomes ||] 

/ [Dz']=4(N-l)(2N-l) [ 1 d8sm 4N - 5 8cos 3 8 (15) 
JQP"- 1 Jo 

The eigenfunctions of the transfer matrix are the Jacobi polynomials pi^-^ 1 ) ( CO s 29) g §. 
Using the integral in the appendix yields 

x? P = ^% ] ^ ) m{1 + 2 > 2N + 2l ' ® (16) 

For N = 2, this is equivalent to the 0(5) magnet. 
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The calculation for a ^ bis more complicated. The reason is that the energy between adjacent 



sites depends now on two angles: the manifold (14) is not a symmetric space. Luckily, the work 



of H allows us to solve this case as well. It is convenient to write coordinates for the quaternionic 
sphere Sp(N)/ Sp(N — 1) in terms of angles 9{ describing an ordinary N — 1-dimensional sphere, 
namely 

q\ = sin On-i ■ ■ ■ sin 02 sin#i 
q2 = sin On-i ■ ■ ■ sin 02 cos 9\ 
q% = sin 9n-i ■ ■ ■ sin 03 cos 62 

qN-l = sin0Ar_i COS0Ar_ 2 

q N = cos 9 N -i 

Then we have 

z 2 i-i = qie 1 ^ 21 - 1 cosoji 
Z2i = qie tip2t sin cjj 

The angles take values in 

< Oi < vr/2 < uji < vr/2 < (fi < 2vr. 

Setting Zi = (0,0, ... ,0, 1) corresponds to all coordinates 9i = ifi = uJi = except ujn = tt/2. 
Then the energy ( |T3| ) is 

( ■ 2 1 1 1 2 1 \ 2 at 
— (asm cj n + ocos tv N )cos a N _ 1 

This indeed depends on two angles, except in the SU (2) gauge-invariant case a = b. In terms of 
these two variables, the measure is 

/ [Dz] =8(N- 1)(2N - 1) / d8sm 4N - 5 8cos 3 6 dwsinwcosw (17) 
Jrs Jo Jo 

where we have dropped the now-unnecessary subscripts and primes. 

The eigenvectors of the corresponding transfer matrix are now labelled by two indices, < 
/' < I < 00. They are written in terms of Jacobi polynomials as || 

cos'"'' 9 (cos 2uj)P( 2N ~ 3 > 1 - 1 ' + V (cos 29) 

One can easily check that they are indeed orthogonal with respect to the measure (|l7|). The 
eigenvalues \u> of the transfer matrix in the Read-Sachdev case a = 0, b = 1 are therefore given 
by 

»7r/2 



A^PS? ) (l)p( 27V - 3 ^' +1 )(-l) = 8(7V-l)(2iV-l) f d0sin 4iv - 5 0cos 3 

Jo 

x dui sin u cos eve 13 cos " " cos " e P^f (cos 2u)P^ N ^ l - l ' +1 ^ (cos 29) 
Jo 

Using the integral in the appendix twice, we find that 



\% S = (-f3) l - l 'T(2N) PAj (1,1') 

j=0 



S 



where 

n _ rg - v + 1 + j)r(|z - jr + 2 + j)r(i / - + i + j) 
jl ' J i!r(2/-2/' + 2 + i)r(i/-|/' + i + i)r(2iv + iz + i/' + i) 

This series be rewritten in terms of a confluent hyper geometric function if desired. For the ground 
state I = V = 0, it simplifies to 

X$ = M(1,2N,0). 

The lowest eigenvalue is identical to that for the SU(2N) magnet, but the general eigenvalues 
are not the same. 



3 Phase transitions as N — > oc 



In this section we study systems with the number of sites L — * oo. This means that the free 
energy follows from the largest eigenvalue Ao of the transfer matrix. We saw in the last section 
that in all cases, Ao is analytic function of the inverse temperature (3 for any finite N. Thus the 
only possibility for a phase transition is if this function develops a singularity as N — > oo. 



3.1 No transition in the 0(N) magnet 

It is useful to study the O(N) magnets first. We need to use the asymptotic formula for Bessel 
functions, valid when 7 is large |8[: 

1 ( fi{\frW-i;- 1 ) 
Iliriy) ~ 7^ (l + y 2 ) 1 ^ 

where 

r 1 = in(i + V1 + 2/ 2 ) - Hv) 

This formula implies that to have a non-trivial large TV limit, we need to take the inverse tem- 
perature (3 to 00 as well, leaving the model in terms of the new variable y = 2(3 /{N — 2). Using 
this formula and some algebra gives the eigenvalue ratios to be 

When we write quantities as a function of y, we mean that the expression is valid up to terms of 
order 1/N. The two-point function is 



(si ■ s i+R ) oc 



R 



so £ is indeed the correlation length. At zero temperature (y — > 00), ^ diverges, just like the 
one-dimensional Ising model. This is the usual zero-temperature behavior in one-dimensional 
classical models. The correlation length does not diverge for any other value of temperature, 
but note that at infinite temperature, £ — » 0. This is a state where the system is completely 
disordered: every site is essentially independent of any other site, because all configurations have 
the same weight in the partition sum. It is also useful to compute the internal energy. This is 
defined as 

d 

U = - — hiA . 
ay 
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Using the asymptotic formula for the Bessel function gives 



N y 

In the N — > oo limit, the internal energy is proportional to N, so it is the energy per component 
which remains finite. 

3.2 Transitions in the SU(N) and Sp(N) magnets 

The SU(N) and Sp(N) magnets each have a phase transition when N — * oo. The lowest eigen- 
value is given by for example 

X^ U{N) = M{l,N,p). 
As detailed in the last section, for the Read-Sachdev case 

\§s = \% U ( 2N \ To find the large N 
behavior, it is useful to examine the differential equation ( |28|) for M(a, N, Ny) directly. Then 
one can see how one can neglect various terms in the equation in various regimes, and then easily 
solve the equation. This sort of analysis is called boundary-layer theory, and the techniques are 



discussed at length in ||. Rewriting (28) in terms of y = /3/N, one has 

y d 2 M x dM 

%w +{1 - v) ik- aM=0 - (18) 

In the large N limit, we can neglect the first term as long as y < 1. This shows that 

lim M(a, N, Ny) = 1 for y < 1. (19) 

N^oo (1 — y) a 



This formula can be verified by using Stirling's formula in the series expression ( |27| ) for M. The 
lowest eigenvalue at large N is therefore 

Af (iV) (y) ~ j^— forj/<l. (20) 
1-y 

Two important results are apparent from this formula. First, there is a singularity at y = 1. 
Second, U su ( N \y) for y < 1 does not grow with N, and so the internal energy per component 
vanishes. A vanishing internal energy is characteristic of infinite temperature. The remark- 
able characteristic of the SU(N) magnet is that this behavior persists all the way from infinite 
temperature y = to a finite temperature y = 1. 

In fact, the correlation length vanishes for all y < 1. This follows from the ratio of the first 
two eigenvalues: 

Af (jV) (j/) f 1 M(2,N + 2,Ny) 

X S U(N \y) [ y) N(N + l) M(l,N,Ny) 

y 



N(l-y) 



for y < 1. 



As N — » oo, this vanishes if y < 1. Most scale-invariant critical points have diverging corre- 
lation length. At the phase transition at y = 1, the correlation length vanishes. We call a 
phase with vanishing correlation length seriously disordered. This distinguishes it from the y > 1 
phase, which is a conventional disordered phase with a finite, non-zero correlation length. In the 
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seriously-disordered phase, the eigenfunction is Pq = 1. This eigenfunction gives equal probabil- 
ities to all configurations. Basically, what has happened is that the energy has been swamped 
by the entropy. For large N, there are so many possible spin configurations that if y < 1, the 
energy term is too small to make a difference. This is apparent in the integral ([ll]). The measure 
favors configurations where the angle 6 between nearest-neighbor spins is near ir/2. The reason 
is as stated before: for a fixed z, there is only one value of z' where z* ■ z' = 1 but there are 
many with z* ■ z' = 0. The energy favors aligned spins, and at y = 1 the energy term is strong 
enough to cause a transition to a phase with finite correlation length. The phase for y > 1 is 
still a disordered phase, but a conventional one, as occurs in the 0(N) magnets. Note that this 
phase transition occurs in the ferromagnetic phase (/? > 0). If (3 < 0, the energy and the entropy 
both favor disorder, so the system is always seriously disordered. 

We have asserted but not yet proven that for y > 1, the internal energy is non-zero and the 
correlation length is finite. The expression (filJ) for M in the large- ./V limit is not valid for y > 1, 



because the first term in the differential equation (18) can no longer be neglected when (y — 1 



is of order l/y/N. To understand y > 1, it is useful to derive a differential equation for 



U 

directly. One has 



%M(l,N,Ny) 
M(l,N,y) 



d M" (M' 

V Ty U = - y W + y {-M 



= UN{y-l)-N + yU z (21) 

where we utilized (|l8|). The internal energy for y < 1 is a solution of this equation with U finite, 
found by neglecting the left-hand side and the last term on the right-hand side. It is 

E^W(y) = -J_ y<l 

y - 1 



The solution of the differential equation (^Tj) for N large and y > 1 comes by assuming U/N is 
finite. Then one neglects the left-hand-side and the second term on the right-hand side, yielding 

USU(N) {y) x _ y 

N y V 

One can find the 1/N corrections to these expressions systematically. For example, one can show 
that the lowest eigenvalue Ao for y > 1 is to next order in N: 

X S o UiN \y)=A N e N ^V~ N y>l (22) 

where An = V^ttN. 

The ratio of the first two eigenvalues for y > 1 follows from a similar computation, yielding 

xf U{N) (y) y M(2,N + 2,Ny) 

x SU(N) {y) N + 1 M (l,N,Ny) 

~ y > 1 

y 



n 



Thus we see that the correlation length indeed goes to zero as y — > 1 from above. This is a 
second-order phase transition: the energy is continuous but its derivative is discontinuous. For 
the Read-Sachdev Sp(N) magnets, the formula is similar: 

x ? s (y) - (vEl „>i 



3.3 Near the phase transition 

The expression for the eigenvalues in terms of the Kummer function M is valid for any value of 
y and N, but the expressions derived for the large- N limit break down when \y — 1| is of order 
(vN). All the terms in the differential equations ( p^j2l| ) need to be included in this region. 
To understand the transition region, we study the physics in terms of the variable 



'N{y-1). 

The differential equation for the lowest eigenvalue (|i~8|) becomes (in the SU(N) case) 



x \ d 2 \ d\ . . 

1 + 7w)l^- x l^- Xo = - (23) 



For N large, we can neglect the term with the 1/yN in it. Thus the subsequent analysis will 
be valid up to terms of order 1/yN, as opposed to equations in the last section, which have 
corrections of order 1/N. Solving this differential equation by plugging in the series 



Ao = E 



3 ' 



j=0 

requires that 

1 

3+ 3+2 3 

Plugging this in and summing the series, the solution near y — 1 for large N is therefore 

A = m e x2/2 + m 1 xM(l,3/2,x 2 /2) (24) 

where the M is indeed our friend the Kummer function. Thus we have shown that the eigenvalue 
Ao = M(l, N, Ny) near y— 1 at large iV is related to a Kummer function with different arguments. 
To fix the values of rriQ and mi , we need to match this onto the value of Ao in the high-temperature 
phase (p0|) and in the low-temperature phase (^2]). From Q, we have the asymptotic formula 

M(l,3/2,z)«iJV z large 



Matching while neglecting terms of order 1/yN yields mo = miy/ir/2. Matching this with the 
numerical result for An given in the last section gives 



m ° = V ~2~ mi = 

The two terms in ( |24| ) cancel when x is large and negative, but add when x is large and positive. 
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With a little work, one can now derive the specific heat 



" d(l/y) 

everywhere, including at y = 1. We find C = 1 for y < 1, C = for y < 1, and C = 1 — (2/7r) for 
y = 1. The traditional critical exponents can be defined when N = oo. For the specific heat, we 
have a = 0, while the correlation length goes to zero logarithmically, so v = as well. 



4 First-order transition for three sites 

In section 2 we found all the eigenvalues and their degeneracies for the SU (N) magnet. This 
means that we can also study phase transitions where the number of sites L is finite as well as 
the infinite-L case studied in section 3. The phase transition persists all the way down to two 
sites. A remarkable result of the mean-field-theory analysis of ]!]] is that for three sites (and only 
for three sites), the transition becomes first-order. We verify this result in this section. 

The partition function for the SU(N) magnet with periodic boundary conditions for any 
number of sites L is given by 



oo 

£ 

z=o 



A(Aj 



where the eigenvalues A; and their degeneracies D\ are given by (12) and ( |10|) respectively. In this 
section we study only the SU(N) magnet, so we will omit the SU(N) superscripts. If one takes 
L — > oo before taking N — > oo, the partition function is dominated by the largest eigenvalue Ao no 
matter what the degeneracy D[ is. However, for finite L, the fact that D[ grows quite quickly with 
I and ./V means that the largest eigenvalue does not necessarily dominate the partition function, 
but instead 

for some Iq which may not be zero. 

All the eigenvalues with finite I (so that l/N ^ as N — > oo) have the same kind of singularity 
as y — > 1 . Thus if the effect of including the degeneracy merely shifts Iq to some finite value, the 
second-order transition remains. To find a first-order transition, we need to study the eigenvalues 
when N — > oo with l/N remaining finite. It is convenient to fix the inverse temperature y and 
study the problem as a function of r, defined as 

r =4- 

At some value r = ro, D(r)\(r, y) L is a maximum; at this value the free energy is minimized (the 
path integral has a saddle point). The behavior of the eigenvalues at non-zero r can be found by 
doing a saddle-point approximation to the integral representation of the Kummer function, or by 
deriving a differential equation for M'/M in the manner of (|2l|). To leading order in N, 

2 

— ln(A(r,y)) = A + y - 1 + ln(2) + r ln(y) + (r + 1) ln(r) 

- ln(y - 1 + A) - (r + 1) In ((r + l) 2 - y + (r + l)A) (25) 

A= ^(r + 1)2 + (y- 1)2-1. 



where 
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At large N and I, the degeneracy behaves as 
1 



— hx(D(r)) = (r + 2) ln(r + 2) - r ln(r) - 2 ln(2) 



(26) 



The degeneracy is a group-theoretical factor and is independent of the inverse temperature y. 
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Figure 1: Plots of the free energy per site at and near the first-order transition for three sites 

A first-order transition occurs if ro jumps discontinuously as y is varied. What can happen is 
that D(r)X(r,y) L can develop another maximum as a function of r. At some value y = yo, one 
can have two peaks at the same place: 

D(n)\{ri,yo) L = D(r 2 )X(r 2 ,y ) L 

i [D(r)\(r,y )\ =ri = ± [D(r)A(r, yo ) L ] r= „ 2 = 0. 

If this happens, the free energy has two minima as a function of r. This is the mark of a first- 
order transition: the internal energy will not be continuous in y because the value of ro jumps 
from r\ to r 2 as y is varied. We find that for L = 3, this indeed happens. In Fig. [I], we plot the 
free energy per site F/3 = — In [D(r)X(r, ?/) 3 ]/(3y) for values of y just below, above, and at the 
transitional value yo = .900312694913. ... At the phase transition, the value of ro jumps from 
to .0197950036795.... This is small, but the transition is definitely first order, as shown in 
PJ]. The transition occurs for y < 1, so the internal energy jumps from to a finite value. No 
transition occurs at y = 1, because the saddle point is already away from r = 0. 

For any other value of L, the transition is second-order as before; no second minimum seems 
to occur. Due to the unwieldiness of the expression (|25|), we have not been able to prove this in 
general, but it is easy to see by looking at the curves numerically. Since the first-order transition 
is so weak for three sites, it would indeed be surprising if the effect of the degeneracies were to 
overcome the behavior of the eigenvalues for larger L. To see this in more detail, let us examine 
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some limits. At large r, the eigenvalue goes to zero as (yr)~ Nr e Nr , while the degeneracy only 
grows as r 2N . Thus ro must be finite. The behavior for small r crucially depends on whether 
y is greater or less than 1. For y < 1, D(0)X(y, 0) = 1, and the slope at r = is negative. 
For y > 1, L>(0)A(y,0) L depends on y and the slope at r = is positive. The simplest behavior 
for y < 1 consistent with these limits is to have D(r)X(y,r) L fall off monotonically from 1 to 
zero. The simplest behavior for y > 1 is for D(r)X(y, r) to rise up to a single maximum at some 
value r = tq, and then fall off to zero for large r. By studying plots of the free energy, it seems 
that this simple scenario is realized for all L larger than 3. For L > 3 and y < 1, ro = 0, so 
the partition function is dominated by an eigenvalue with finite I. As y is increased past 1, ro 
becomes non-zero. There is only a single peak, so ro varies continuously with y, so the minimum 
value of the free energy is also varies continuously. This means that for L > 3 the only phase 
transition is a second-order one at y = 1. 

We are very grateful to Shivaji Sondhi for many interesting conversations and for collaboration 
on We are also grateful to R. Moessner and N. Read for helpful conversations. The work 
of P.F. is supported by NSF Grant DMR-0104799, a DOE OJI Award, and a Sloan Foundation 
Fellowship. O.T. is supported by NSF Grant DMR-9978074. 

A Jacobi polynomials and confluent hyper geometric functions 

The Jacobi polynomials P^ ,u '{x) can be expressed in terms of a Rodrigues formula 

(i _ + x y P ^\ x) = ^(A}j l [ ( i _ x y+^! + x y+» 

They are defined to be orthogonal with respect to the measure dx{l — + x) u : 

dx(l - xf(l+x) u P l M (x)P^' v \x) =0 l^k 
Many results on Jacobi polynomials can be found in |10[ |. One integral we need is (l 

J -i U T(a - v — I + l)T(fi + a + I + 2) 



Be aware that there are (different) typos in this formula in both []TT| and in |T2]. Another useful 
result is 

r ( u ,u )(1) r(z + M + i) 
lJ r(z + i)r(/x + 1) 

Using these two relations we evaluate the integral 



dx (1 - x)"(l + ^e^+^p/^Or) 
i=o J ' 

= V /? .7 2 M+^+:/+l + j + 1 ) r (/ i + j + 1 

2-^ P l\r(n - 1 -U Firr,< 4- ; y 4- i 4- 



i!r0'-z + i)r(// + i/ + j + z + 2) 
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P^ v \l) Tip + 1) l f^{2(i) k 2^ +l+1 r ^ + k + l + 1 ) 



k\T(fi + u + k + 2l + 2) 

= P^\l) (3 l 2^ +i r ^ ) ^+| + 2) 1) Af(, + 1 + 1,^ + 21 + 2,20) 

where k = j — I and we use the fact that 1/T(j — I + 1) = for j — I a negative integer. The 
function M(a, b, z) is known as Kummer's series: 

r (°) ^ kl F(b + k) 

The Kummer function is a confluent hypergeometric function because it can be obtained by 
taking a limit of the hypergeometric function 2 -Pi where two singularities coincide. It satisfies 
the differential equation 

d 2 M s dM 
z— r + (b-z)—-aM = (28) 
az z az 

M is an analytic function of a, b and z; the only way to get singularities is to take some or all 
of these parameters to infinity. A fact useful for O(N) magnets is that the Bessel function I 7 (z) 
can be written as 

W = 2 "r(7 + i) e ~* M(7 + 1/2 ' 27 + 2z) (29) 
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